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INCLUDING THE EFFECTS OF GAMMA 
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SUMMARY 

A method for calculating two-dimensional inviscid nozzle contours for a real gas or 
an ideal gas by the method of characteristics is described. The method consists of a mod- 
ification of an existing nozzle computer program. The ideal-gas nozzle contour can be 
calculated for any constant value of gamma. Two new methods of calculating the center- 
line boundary values of the Mach number in the throat region are also presented. The use 
of these three methods of calculating the center-line Mach number distribution in the throat 
region can change the distance from the throat to the inflection point by a factor of 2.5. 

The effect of ratios of specific heat (gamma) equal to 1.6667, 1.4, 1.25, and 1,1 on 
two-dimensional nozzle contours with constant area ratios and nearly constant lengths 
showed that a decrease in gamma resulted in a much more gradual opening of the nozzle 
contour from the throat into the expansion region. For similar area ratio and length con- 
straints and for the same values of gamma, calculations of axisymmetric nozzle contours 
showed a trend similar to that of the two-dimensional nozzle contours; however, the effect 
of gamma on the axisymmetric contours was not nearly as pronounced as the effect on the 
two-dimensional contours. The effect of gammas of 1.6667, 1.4, 1.25, and 1.1 was inves- 
tigated for two-dimensional and axisymmetric inviscid nozzle contours with a constant 
exit Mach number. The effect of decreasing the gamma in the calculation of the two- 
dimensional nozzle contour resulted in a more gradual opening of the nozzle from the 
throat into the expansion region. The effect on axisymmetric contours for calculations 
made with the same four values of gamma was similar to the effect found for the two- 
dimensional contours; however, the effect on the axisymmetric contours was not as pro- 
nounced as the effect on the two-dimensional contours. For the two-dimensional nozzle, 
contours with a constant exit Mach number, the ratio of the exit height to throat height 
increased sharply with a decrease in the value of gamma. For a fixed throat size, the 
length of the two-dimensional nozzle decreased significantly with an increase in the value 
of gamma, 

A user's guide is presented for input to the computer program for both the two- 
dimensional and axisymmetric nozzle contours. 



INTRODUCTION 


The calculation of supersonic and hypersonic nozzle contours by the method of char- 
acteristics has been used in the design of turbines, rockets, and wind tunnels. A review 
of many of the methods of supersonic nozzle calculation is given in reference 1 , while 
additional methods are presented in references 2 to 9. This paper presents some of the 
major modifications made to an inviscid nozzle computer program since it was presented 
in reference 10. The modifications include a new technique of obtaining the starting center 
line and a method for calculating two-dimensional nozzle contours for a real gas or an 
ideal gas by using a modified version of the axisymmetric program presented in refer- 
ence 10. In addition, a user’s guide is given in order that the proper input can be used 
to select one of the 24 options that may be used in the calculation of an inviscid nozzle 
contour. Other input parameters used to form boundary conditions for the characteristic 
mesh are presented in the form of curves for ratios of specific heat (gamma) of 1.6667, 

1.4, 1.25, and 1.1. Sample nozzle contours are also presented for the above four values 
of gamma and for parts of the nozzles which use a high temperature combustion gas. Por- 
tions of the nozzle contours calculated using the high temperature combustion products 
are presented for the maximum turning angles (0 C ) of 7.5°, 12°, 15°, 18°, and 20° and 
for the throat regions with three different center-line Mach number distributions. 

SYMBOLS 

A,B,C,D,E coordinate points for starting nozzle calculations (cf. sketches (a) and (b)) 

A coefficient defined by equation (1 3) 

A'jB’jC’ coefficients defined in equation (15) 
a velocity of sound 

C,D coefficients defined in equation (1) 

* ~\ 

Ify stagnation enthalpy 

L nondimensionalizing scale factor for nozzle contours 

M Mach number 

p pressure 

r distance measured from the source point in a radial flow field 


2 



r cr 

s 

u 


u l 

V 

W 


X, 


me 


x 


y 

y 

e 

d i 

®c 

P 

P 

P 

P 


distance measured from the source point to a sonic sphere in radial flow 

distance along a Mach line 

velocity component in the x direction 

limiting velocity 

velocity component in the y direction 
limiting velocity ratio, u/u, 

incremental change in x distance defined by equation (14) 

distance along nozzle plane of symmetry (or axis) with x = 0 at 
nozzle throat 

distance measured perpendicular to nozzle plane of symmetry (or axis) 
ratio of specific heats 
flow angle, rad 

integrated values of flow angle as given by equation (A6) 
maximum turning angle, deg 
Mach angle, sin"^^ 
density 

stream function 

nondimensional stream function 


Subscripts: 

A,B,C,D,E properties evaluated at points defined by sketches (a) and (b), with A 
referring to point evaluated at the nozzle throat 
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lim value of the stream function that defines the inviscid nozzle contour 

t total conditions in the stagnation chamber 

°° free-stream uniform flow conditions 

METHOD OF CALCULATION 

Determination of Center-Line Mach Number Distribution 

A detailed method of calculating axisymmetric inviscid nozzle contours was pre- 
sented in reference 10. This method included a means of obtaining the boundary condi- 
tions for the construction of a characteristic mesh which used a linear Mach number dis- 
tribution in the throat region. The method for obtaining the boundary conditions presented 
in this report is a modified version of the method presented in reference 10 and can be 
applied to the calculation of either two-dimensional or axisymmetric inviscid nozzle con- 
tours. In order to construct a characteristic mesh by the method of reference 10, the 
Mach number, the flow angle, and the value of the stream function must be known along 
the center line in the expansion region, along the radial-flow region boundaries, and along 
the final Mach line. The new method presented in this report starts the characteristic 
mesh at the final Mach line and calculates along Mach lines which start at the nozzle cen- 
ter line. The calculation ends at a point in the throat region where the Mach number is 
slightly greater than 1.0. On the center line the stream function and flow angle are zero 
and along the final Mach line the Mach number is constant, the flow angle is zero, and 
the stream function is easily calculated along the final Mach line as a function of the 
y coordinate. Thus, the only boundary condition to be determined is the Mach number 
distribution along the center line. 

As presented in reference 10, the method of calculating the initial flow conditions 
for starting the characteristic mesh divided the inviscid nozzle flow into four regions as 
shown in sketch (a). The flow in region TV is a uniform parallel flow at the design Mach 
number of the nozzle. The triangular region HI (BCD) is a radial-flow region which 
defines the flow along the boundaries B to C and C to D. In reference 10 a linear Mach 
number distribution is assumed from point A to B, where the Mach number is defined by 
the radial-flow region. The boundary from D to E is the final Mach line and has a con- 
stant Mach number and a 0° flow angle. Regions I and n use the flow properties on the 
boundaries A to B to C and C to D to E, respectively, to start the construction of the char- 
acteristic mesh. The calculation starts at point D and moves up the Mach line that is just 
upstream of the final Mach line (Mach lines of the first family shown by dashed lines in 
sketch (a)) until the integrated value of the stream function (along the Mach line) is reached. 
It can be seen from sketch (a) that as the calculation proceeds from point D to point C the 
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Mach lines (dash lines) decrease in length until just before point C is reached the Mach 
line is extremely short. After the Mach line closest to point C is constructed, the cal- 
culation returns to a Mach line just upstream of point B in region I. Sketch (a) indicates 
that the last Mach line calculated in region II is much shorter than the first Mach line 
calculated in region I. 

These two different length Mach lines are used to calculate adjacent inviscid bound- 
ary (contour) points which are determined by a numerical integration of the stream func- 
tion along the Mach lines. The numerical integration of the stream function along two 
different length Mach lines results in a discontinuity in the inviscid contour as demon- 
strated at the inflection point (point C) in figure 1. The discontinuity at point C is believed 
to result from a cumulative error in the numerical integration of the stream function, cou- 
pled with a small cumulative error inherent with a numerical finite -difference calculation 
with successive approximations. From figure 1 it should be noted that the error in the 
y ordinate caused by the discontinuity is less than 2 percent; however, this discontinuity 
can be removed by using a different method of obtaining the boundary condition for the 
characteristic mesh. 

The discontinuity at the inflection point can be removed by starting all the mesh 
points (Mach lines) from the center line (points D to B to A) as shown in sketch (b). The 
center-line Mach number distribution from point D to point B is calculated using the same 
radial-flow equations used previously in sketch (a) along B to C and C to D. The center- 
line Mach number distribution is obtained by scribing an arc with a radius of from 

point P ? down to point P on the center line (see sketches (a) and (b)) where the Mach num- 
ber at point P' equals that at point P. Thus, the entire characteristic mesh can be started 
from the center line with Mach lines gradually decreasing in length, and the discontinuity 
at the inflection point found with the method of reference 10 can be avoided. 
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In reference 10 the method of determining the Mach number distribution from point A 
(nozzle minimum) to point B used a linear Mach number distribution. This linear distri- 
bution of Mach number from A to B produced a finite Mach number gradient with respect 
to axial distance at the throat. It was pointed out in reference 1 that in order to have a 
straight sonic line normal to the axis at the nozzle minimum, the Mach number gradient 
must be zero at the sonic point. It was also pointed out in reference 1 that experimental 
results have shown that the best flow is obtained when small maximum turning angles are 
used in conjunction with the requirement that ■ - = 0 at the throat; this requirement 

1v) 

results in a relatively long smooth contour that approaches zero curvature at the sonic 
point. Therefore, in order to improve the quality of the flow in the throat region the 
method of obtaining a Mach number distribution along the center line from point A to B 
was changed from a linear expression (ref. 10) to 
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(Referred to as ABKEY = 1) 
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The boundary conditions at 
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Equations (2) and (3) cannot be solved until a value for is determined. 

Between the two end points of the throat region i > it is required that 
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The limi tation given by equation (8) yields an increasing Mach number as increases. 

The condition of equation (9) gives a Mach number gradient which increases monotonically 
from a value of zero at the throat to a value of f— at B» 
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Equations (10) and (11) indicated that in order for the conditions of equations (8) 
and (9) to be satisfied, that C and D must be greater than zero. Positive values of 
C and D are assured if a parameter found from the terms in equations (2) and (3) has 
the following limits 



Equation (12) implies that the value of 

V r cr/ B 


is found from 



( 12 ) 


(13) 



Two methods of calculating a center -line Mach number distribution from A to B 
(see sketch (a)) have been described. The first method stipulated a linear Mach number 
distribution from A to B with a finite slope of /v y at = 0. The second method 

F©J b cr 

used a third order polynomial (eq. (1)) to describe the distribution of Mach number as a 
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function of The boundary conditions stipulated that r = 0 at = 0. The 

rcr d fe) 

use of the second method resulted in a throat region that was two to three times longer 
than the region calculated from the linear variation of Mach number. 

If the nozzle design requires a shorter throat region from A to B than can be calcu- 
lated using a linear Mach number distribution, a third method for finding the Mach num- 
ber distribution from A to B can be used. The third method used a second order poly- 
nomial given by 

M = A* + B + C \r^f (Referred to as ABKEY = -1) (15) 

where 


A’ = 1.0 (16) 



The boundary conditions are: 
At 
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The stipulation in the use of equation (16) is that any finite value of 
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) that is smaller than the value found from the previously discussed 
V r cr/ B 

method and also smaller in value than the method of reference 10 (i.e., linear Mach num- 
ber distribution). When the slopes at points A and B are equal, equation (19) reduces to 


the linear Mach number distribution of reference 10. When 


times greater than 


the 


(*)» 


dM 


dl r~ 

L cr 


dM 


cr, 


is two or three 


, the is two-thirds or one-half, respectively, of 


cr/ 


B 


J B 


distance found using the linear Mach number distribution from A to B. 


cr, 


After the value of MM is determined, equations (15) to (18) are solved for the 
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Mach number distribution along the center line from A to B. The values of between 
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point B and D (including point D) are changed to values by adding X inc to the 
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RESULTS AND DISCUSSION 


Two-Dimensional Nozzle Calculations 

The design of a two-dimensional real -gas nozzle contour was required as part of 
the design of a thermal protection system test facility at the Langley Research Center; 
this facility burns a mixture of air, oxygen, and methane in a settling chamber and expands 
the products of combustion in the nozzle. These products of combustion have an effective 
value of gamma of approximately 1.25. The method used for calculating two-dimensional 
nozzle contours is given in the appendix. The combustion facility will be used to test 
space shuttle panels under the conditions of heating expected in flight. The mechanical 
design criteria necessitated that the nozzle have as short a length as possible (to avoid 
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excessive heat transfer) within the limits of the two-dimensional method of calculation 
described in this report. In order to assess the effect of the maximum turning angle on 
the length of the nozzle, real- gas inviscid nozzle contours were calculated for maximum 
turning angles of 7.5°, 12°, 15°, 18°, and 20° and are shown in figure 2, The nozzle 
contours were calculated with a linear Mach number distribution on the center line of 
the expansion region (A to B, sketch (b)) and for a real-gas Mach number of 3.777. The 
Mach number was selected from the real-gas flow properties at a stagnation pressure 
and enthalpy of 10 atm and 5.307 x 10 6 joules/kg in order to give an inviscid area ratio 
of 25. The results in figure 2 showed that as the maximum flow angle increases the 
length of the nozzle decreases. For instance, an increase in the maximum flow angle 
from 7.5° to 15° reduces the length of the nozzle to slightly half the original length. 

The nozzle can also be shortened by selecting the Mach number distribution along 
the center line (A to B, sketch (b)); the ABKEY = -1 gives the shortest expansion region. 
Figure 3 illustrates the effect of using the three different center-line Mach number dis- 
tributions; where ABKEY = 1 is the distribution given by equation (1), ABKEY = 0 is 
the linear distribution used in reference 10, and ABKEY = -1 is the distribution given by 
equation (15). The nozzle calculations in figure 3 were made for a combustion facility 
using real-gas flow properties, with a design Mach number of 3.777. The expansion 
region with = 0 at the throat (ABKEY = 1) has the greatest distance from the throat 

d © 

to the inflection point C ~ 51j . The linear distribution (ABKEY = 0) shortens the 

distance from the throat to the inflection point to almost half the distance calculated by 
the ABKEY = 1 method. The use of the ABKEY = -1 method compared to the 
ABKEY = 1 method reduces the length of the nozzle from the throat to the inflection 
point by a factor of approximately 2,5. Ideally, a nozzle should have a zero Mach number 
gradient at the throat as given by the ABKEY = 1 design. However, the very gradual 
opening of the throat region of the nozzle using the ABKEY = 1, as seen in figure 3, pre- 
sents design and fabrication problems as well as an excessive heat load to the nozzle. 
Problems such as these dictate that the ABKEY = -1 design method should be used for 
a nozzle that operates at high enthalpy stagnation conditions. The combustion products 
have an effective gamma of approximately 1,25 which coupled with the ABKEY = 1 
method of calculating the throat region results in a gradual opening of the throat region as 
previously noted in figure 3. The very gradual opening of the throat region for the 
ABKEY = 1 method would be accentuated in the design of nozzle contours for a gas that 
has a gamma of approximately 1.1 (similar to tetrafluoromethane (CF^)]. Thus, for noz- 
zles with small values of effective gammas (i.e., y ~ 1.1), both the maximum turning angle 
and the type of throat region must be carefully chosen to satisfy certain mechanical and 
heat transfer design criteria. 
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The Effect of Gamma On Nozzle Contours 


The results in figure 3 indicated that for a gamma of approximately 1,25 the nozzle 
had a very gradual opening in the expansion region. In order to study the effect of gamma 
on nozzle contours, two-dimensional calculations were made for gammas of 1.6667, 1.4, 
1.25, and 1.1, these gammas represent the value of gamma for helium, ideal-gas air, 
combustion products (or high temperature air), and tetrafluoromethane, respectively, and 
are shown in figure 4. The calculations were made for a constant value of gamma using 
ideal-gas egressions. The four contours were calculated for a constant throat -to-exit 
area ratio of 25 and a constant maximum turning angle of 12°. In order to get a constant 
area ratio for the four values of gamma, the design Mach number was varied from 7.11 to 
3.49 from the highest to lowest value of gamma, respectively. In addition to the constant 
area ratio, the four contours had the same total length at an x/L value of about 64. For 
two-dimensional nozzles of a given area ratio and given length, the effect of decreasing 
gamma (see fig. 4) results in a more gradual opening of the nozzle from the throat into 
the expansion region. 

A comparison of four axisymmetric nozzle contours calculated for gammas of 1.6667, 
1.4, 1.25, and 1.1 is shown in figure 5. The comparison in figure 5 is for the same val- 
ues of gamma, the same Mach numbers, and the same constant area ratio as was used in 
figure 4. The effect of gamma on the axisymmetric nozzle contours is similar to the 
effect found in figure 4; for example, there is a more gradual opening of the expansion 
region as the value of gamma decreases. However, for the axisymmetric contours in 
figure 5, the variation in the expansion region is not as pronounced as was found for the 
two-dimensional contours in figure 4. 


Two-dimensional nozzle contours calculated at a constant exit Mach number of 5.0 
for values of gamma of 1,6667, 1.4, 1.25, and 1.1 are shown in figure 6, The four nozzle 
contours had a fixed throat size ^ = 0.24^ for the purpose of facilitating a comparison 
of the effect of gamma on the contours. The nozzle contour in the expansion region just 
downstream of the throat begins to open more rapidly as the value of gamma increases 
from 1.1 to 1.6667. The area ratios for the four contours have values of 597.7, 69.22, 
25.00, and 9.804 for values of gamma of 1.1, 1,25, 1.4, and 1.6667, respectively. For 
two-dimensional nozzles, the variation in area ratio is accentuated because the area ratio 
varies directly as the ratio of the exit height to the throat height. Thus, for a fixed throat 
size the nozzle calculated for y = 1.1 is large relative to the nozzle calculated for 
y a 1.6667, The ratio of nozzle lengths for the nozzle designed for y - 1.1 to the one 
designed for y ~ 1.6667 is approximately 80. 


Axisymmetric nozzle contours calculated for values of gamma of 1,6667, 1.4, 1.25, 
and 1.1 and for a constant exit Mach number of 5.0 are shown in figure 7, Again, the four 
contours had the same throat radius 0,24, i.e., the same throat height as that used 
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for the two-dimensional contours in fig. 6j for the purpose of facilitating a comparison of 
the effect of gamma on the contours. Similar to the two-dimensional contours, the con- 
tour in the expansion region begins to open more rapidly as the value of gamma increases 
from 1.1 to 1.6667. The area ratios for the four contours are the same for the respective 
values of gamma as the two-dimensional nozzles in figure 6. However, the variation in 
exit radius for the axisymmetric nozzles is not as great as that for the two-dimensional 
nozzles. For the axisymmetric contours, for a given area ratio, the variation in exit 
radius varies as the square of the radius. Thus, the variation in exit radius is much less 
for the axisymmetric nozzles than it is for the two-dimensional nozzles. The ratio of noz- 
zle lengths for an axisymmetric nozzle designed for y = 1.1 to the one designed for 
y = 1.6667 is a factor of approximately 11 compared to the factor of approximately 80 
found for the two-dimensional nozzle contours. Again, for a fixed throat radius, as was 
found for the two-dimensional contours, the size of the nozzle for y = 1.1 is much 
larger than the nozzle calculated for y = 1.6667. 

CONCLUDING REMARKS 

A method for calculating two-dimensional inviscid nozzle contours for an ideal or 
real gas is described. The method of calculation is similar to a previously reported 
method for calculating axisymmetric nozzle contours. A modified version of a previous 
method for obtaining the boundary conditions for starting the characteristic mesh of both 
the two-dimensional and axisymmetric inviscid contours is presented. In addition, two 
other methods of obtaining the Mach number distribution in the throat expansion region 
are presented. If required, the program can print out the entire flow-field properties as 
determined by a solution of the characteristic mesh. Normally the program just prints 
the inviscid nozzle coordinates with the flow properties (necessary for a boundary-layer 
calculation) given at the edge of an inviscid contour. 

Calculations were made to minimize the length of a two-dimensional real-gas invis- 
cid nozzle contour to be used in the design of a combustion facility that would test space 
shuttle panels. A variation of the maximum turning angle for two-dimensional nozzles 
showed that an increase in the flow angle from 7.5° to 15° reduced the length of the nozzle 
by approximately a factor of 2. The two new methods (designated ABKEY = -1 and 
ABKEY = 1) of determining the center-line Mach number distribution in the throat expan- 
sion region gave different center-line lengths from the throat to the inflection point. One 
of the methods shortens the nozzle (ABKEY = -1) whereas the other method (ABKEY = 1) 
gives a longer expansion region and also gives the proper boundary conditions for a 
straight sonic line perpendicular to the nozzle axis. The ABKEY = -1 (short -throat 
expansion region) method shortens the center-line distance from the throat to the inflec- 
tion point by a factor of 2.5 compared to the ABKEY = 1 method (long-throat expansion). 
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The effect of ratios of specific heats (gamma) of 1.6667, 1.4, 1.25, and 1.1 was 
investigated for two-dimensional and axisymmetric inviscid nozzle contours with variable 
exit Mach numbers, constant area ratios, and constant lengths between the throat and the 
exit. The effect of decreasing the gamma in the calculation of the two-dimensional nozzle 
contours results in a more gradual opening of the nozzle from the throat into the expansion 
region. The effect on axisymmetric contours for calculations made with the same four 
values of gamma was similar to the effect found for the two-dimensional contours; how- 
ever, this effect on the axisymmetric contours was not as pronounced as that on the two- 
dimensional contours. 

Two-dimensional nozzle contours which were calculated for gammas of 1.6667, 1.4, 
1.25, and 1.1 with a constant exit Mach number, showed that a decrease in gamma resulted 
in a much more gradual opening of the nozzle contour in the area of the expansion region. 
For the contours with a constant exit Mach number the ratio of the exit height to throat 
height increased sharply with a decrease in gamma. For a fixed throat size the length 
of the two-dimensional nozzle decreased significantly with an increase in the value of 
gamma. Axisymmetric nozzle contours, calculated with a constant exit Mach number and 
for the same four values of gamma, showed that a decrease in gamma gave a more grad- 
ual opening of the nozzle contour in the expansion region. This trend was similar to the 
two-dimensional nozzle contours; however, for the axisymmetric nozzle contours it was 
not nearly as pronounced. 

A user’s guide is presented in order to input the various options for two-dimensional 
and axisymmetric inviscid contours. 

Langley Research Center 
National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 6, 1975 
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APPENDIX 


ANALYSIS FOR CALCULATING TWO-DIMENSIONAL INVISCID 

NOZZLE CONTOURS 


Characteristic Equations 

The procedure for calculating an axisymmetric inviscid nozzle contour for a real 
gas is described in reference 10. The method of calculation used for the two-dimensional 
(2D) inviscid nozzle contour is very much the same as the method described in refer- 
ence 10. The new technique of obtaining the center-line Mach number distribution and 
the three throat regions, both of which have been previously described in this report, are 
incorporated into both the two-dimensional method described herein and the axisymmetric 
method of calculation described in reference 10. The method of calculation allows for a 
complete description of the flow field as defined by the characteristic mesh and if required 
the entire flow-field properties can be printed by the program. 

The assumption of a steady, irrotational flow results in the basic continuity and 
Euler’s equations given as (see ref. 3, p. 19)t 


3(pu) 9(pv) _ 

ax + ay 


(Al) 


U M + v 9u = _iiE 
ax ay p dx 


dv av 1 8 P 
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From equations (Al) and (A2) the two-dimensional characteristics equations are 
derived in the form (see ref. 3, ch. 2): 
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The characteristic equations are then written in a difference form and solved in a 
stepwise procedure by 3-point successive approximation. Further details are given in 
appendix A of reference 10. 

Flow Properties Along Flow-Region Boundaries 

As previously discussed in this report and in reference 10, the boundary conditions 
along the center line (A to D) and along the final Mach line (D to E) are required to start 
the characteristic mesh. (See sketch (a).) The region on the center line from B to D is 
determined by radial flow. Two-dimensional radial or source flow is defined as a flow in 
which all the flow properties are a function only of the distance r from a fixed point in 
the two-dimensional plane. The general equation for two-dimensional radial flow derived 
from "one -dimens ion gas dynamics" is 

(M 2 -l)^.-^=0 (A5) 

' W r 


For the radial-flow region the equations for the integrated total expansion angle 0j and 

the integrated distance — as taken from equations (A3) and (A5), respectively, are 

r cr 
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Calculation of Streamline Along an Inviscid Boundary 

The fact that the boundary for the characteristic mesh no longer goes from B to C 
and from C to D (see ref. 10) does not mean that point C does not exist. Point C is still 
the inflection point on the inviscid boundary where the flow has its maximum turning angle. 
Also, as a carryover from the method described in reference 10, point C is where the 
radial-flow stream function ^ is evaluated in order to establish the limiting value 
The value of is used to determine the inviscid edge point on the first -family char- 

acteristic. The differential form of the two-dimensional radial-flow equation is 


di// = pur d0 (A8) 

It is convenient to define a nondimens ional stream function as 

V- ' ( A9 > 

Wcr 
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The resulting expression for the radial-flow stream function evaluated at point C is 

r ^ (A10) 


b 


= P.1 

im I p t u z r cr J c 


The general differential form of the stream function for two-dimensional flow is 

dip = pa ds (All) 

where a is the local speed of sound, and ds is the differential distance along a Mach 
line. When equations (A9) and (All) are combined, the nondimensional form of the stream 
function is given by 


\p = & /T ^ £- W sin jx 

Jq p t V cr / 


(A12) 


The value of ip is calculated at each point along the first-family characteristic until 

a value which exceeds is reached. The flow properties and coordinates along the 

edge of inviscid contour are determined by interpolation between the i p values which bound 

^lim* ^ Preliminary design of a two-dimensional or axisymmetric nozzle contour, it 

is convenient to know the integrated value of the flow angle 0 T and the value of —— over 

i r cr 

the boundaries of the radial-flow region. (See sketches (a) and (b).) The values of 0- { 
and as a function of Mach number for gammas of 1.1, 1.25, 1.4, and 1.6667 are shown 


cr 


in figures 8 and 9 for two-dimensional radial flow and in figures 10 and 11 for three- 
dimensional radial flow. The expressions for 0 j and 


for two-dimensional flow are 


cr 


given in equations (A6) and (A7), respectively, and the two expressions for axisymmetric 
flow are given in reference 10. The curves in figure 8 or 10 may be used to obtain an 
estimate of the nozzle length for a given maximum turning angle 0 £ and a given design 


Mach number. The value of 0 


I,B 


and 0 


I,C 


can be found from 


0 I,B “ 0 I,D " 20 C 
0 I,C = 0 I,D " 0 C 

where 0j g is evaluated at the design Mach number. The Mach numbers at point B 
and C can be found with 0j g and 0j q from figure 8 or 10. From the values of the 
Mach number at point B and C an assessment can be made of the Mach number distribu- 
tion in the expansion region. The length of the nozzle can be determined by first evalu- 
ating xpn m at point C from equation (A10), then setting ip of equation (A12) equal to 
^lim anc * solving equation (A12) for at sonic (nozzle -throat) conditions. The value 

of 


cr 


cr 

at the throat is used to obtain a scale factor for the desired throat size. The 


cr 


scale factor is then used to convert ( -±— ) from figure d or 11 to a physical dimension 


D 


which is approximately the distance from the nozzle throat to point D; i.e., the point on 
the nozzle center line where the final Mach line originates. 
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APPENDIX 


$N1 

DOPT 

GOPT 

GAM 

CLINE 

EMTAB 

WTAB 

RHORT 

PPG 

N 

THETC 

EMD 

ABAR 


Program Input Format 
Card 1, starting in column 2. 

2 for two-dimensional contour. 

3 for axisymmetric contour. 

1 for equations for gas properties; i.e., 1 + (~^\, £-= (l - W 2 

y w 2 v ~ 1 \M 2 / p t 

and £-= (l - W 2 ) 7 1 

2 for gas tables of l/w 2 , i/m 2 , pjp^, and p/p^.. 

Value of y if GOPT = 1. 

0 if GOPT = 2. 

1 if center-line Mach number distribution input on $N2 cards. 

0 if radial-flow equations are used to determine center-line Mach number 
distribution. 

Table of i/m 2 , increasing values, limit of 200. Omit if GOPT = 1. 

Table of l/w 2 , increasing values, limit of 200. Omit if GOPT = 1. 

Table of P/P^, increasing values, limit of 200. Omit if GOPT = 1. 

Table of p/p t , increasing values, limit of 200. Omit if GOPT = 1. 

Number of values in a gas table; i.e., number of values in i/m 2 table. 
Omit if GOPT = 1. 

Maximum turning angle at point C, in radians. Omit if CLINE = 1. 
Design Mach number, value at point D. Omit if CLINE = 1. 

Value of A used with ABKEY = 1 for calculation (eq. (13)). 

\ r» y » I 
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ABKEY -1 for shortest throat expansion region (requires value of DMDXA). 

0 for linear center- line Mach number distribution in the throat expansion 
region. 

1 for longest throat expansion region (third order polynomial for Mach number 
distribution of the center line). 

Omit if CLINE = 1. 


DMDXA 


Value of 



at point A (throat). If ABKEY = -1, omit if CLINE = 1. 


DELW The integration interval (dW) for the radial-flow equations (eqs. (A6) and (A7) 
and ref. 10). Omit if CLINE = 1. 

RLIM The value of r/r cr above which the r / r CT spacing (Ar/r cr j is kept within 
a limit defined by DRORCR. (This must be determined from an initial trial 
calculation of the radial -flow region.) Omit if CLINE = 1. 


DRORCR 


Limit on the r/r cr 
CLINE = 1. 


spacing 



after RLIM is reached. 


Omit if 


RERR 

DDW 

DXDE 

NABB 

MAX 

$ 


Error limit on the value of DRORCR. Omit if CLINE = 1. 

Change in DELW required to keep DRORCR within the error limit RERR. 

Omit if CLINE = 1. 

Delta along DE. Omit if CLINE = 1. 

r cr 

Number of points from A to B, Omit if CLINE =1. 

The maximum number of boundary points on the center line and the final Mach 
line, set at 1500. Omit if CLINE = 1. 

Start in column 2. 


For CLINE = 1: 


$N2 Start in column 2. 


XC value on the center line at the origin of the final Mach line. Omit if 

r cr 

CLINE = 0. 
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EMC Mach number along the final Mach line originating at XC. Omit if CLINE = 0. 

DXBC The -2- increment along the final Mach line. Omit if CLINE = 0. 

r cr 

PSID Value of stream function used to define inviscid contour. Omit if CLINE = 0. 

$ In column 2. 

For the CLINE = 1 option each point on the center line from the point just 
upstream of the final Mach line to throat must have a value of XC and EMC on a card with 
the following format: 

$N2 Column 2. 

XC Point on the center line. 

EMC Mach number corresponding to XC. 

$ Column 2. 

After the center-line Mach number distribution has been input a final card is input: 

$N2 Column 2. 

XC A value of 0. 

EMC A value of 0. 

$ Column 2. 

The next set of inputs are: 

$N1 Column 2. 

NAPRX Set equal to 3, It is the number of approximations used in the solution of the 
characteristic mesh points. 

PRINT 1 print characteristic network, 

0 print inviscid wall points only. 
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AREA 1 starts case. 

NN Number of values in table of l/M^. Omit if GOPT 

$ Column 2. 

The final card is: 

$N1 Column 2. 

AREA 2 ends case. 

$ Column 2. 
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(a) End of expansion region, including inflection point and start of turning region. 

Figure 1,- A comparison of two-region method of calculation with the one region of 
method of calculation for an inviscid nozzle contour. 
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Dimensionless distance from nozzle throat, x/L 


Figure 2.- The effect of various maximum turning angles on two-dimensional real-gas nozzle contours 
(ABKEY = 0), (p t = 10 atm, H t = 5.307 x 10 6 joules/kg), and (y ~ 1.25). 
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Figure 3.- The effect of three different center- line Mach number distributions in the expansion region of a real- gas 

two-dimensional nozzle contour (y * 1.25). 
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Figure 10.- The integrated flow angle as a function of Mach number for axisymmetric flow and various ratios 

of specific heats. 
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